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Abstract 



Based on the short-time dynamic scaling form, a novel dynamic approach 
is proposed to tackle numerically the Kosterlitz-Thouless phase transition. 
Taking the two-dimensional XY model as an example, the exponential diver- 
gence of the spatial correlation length, the transition temperature Tkt and 
all critical exponents are computed. Compared with Monte Carlo simulations 
in equilibrium, we obtain data at temperatures nearer to Tkt- 
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The Kosterlitz-Thouless (KT) phase transition is an important kind of phase transitions 
in nature |l|J^]. When the temperature approaches the transition temperature Tkt from 
above, the spatial correlation length diverges exponentially, rather than by a power law in 
a second order phase transition. Below Tkt, the system remains critical in the sense that 
the spatial correlation length is divergent. No real long range order emerges in the whole 
temperature regime. Important examples of systems with a KT transition are the classical 
XY-type models, quantum Heisenberg models, hard disk models and other relevant fluid 
systems as well as field theories. 

It is well known that due to the exponential divergence at the transition temperature, 
numerical simulations of critical systems with a KT transition suffer severely from critical 
slowing down. For example, to compute the spatial correlation length of the two-dimensional 
classical XY model, even with the cluster algorithm and the over-relaxed algorithm one has 
only reached the temperature T = 0.98, which is still fairly far from Tkt estimated to be 
around 0.89 to 0.90 If some quenched randomness is added to the system, e.g. in the 

fully frustrated XY model, the situation becomes even more complicated |6|-[T0 . 



Recently much progress has been made in critical dynamics. It is discovered that univer- 
sal dynamic scaling behavior emerges already in the macroscopic short-time regime, after a 
microscopic time scale t mic [TT|-|T9"| . More interesting and important is that the static expo- 
nents originally defined in equilibrium enter the short-time dynamic scaling. This provides 
a possible way for extracting these exponents from the short-time dynamic scaling behavior 
||19|| . Since the measurements are carried out in the macroscopic short-time regime, the 
method is free of critical slowing down. 

Such a short-time dynamic approach has been systematically investigated for critical 
dynamic systems with a second order phase transition. It has been first verified in the 
simple Ising and Potts model |2(]f2Ijand recently applied successfully to general and complex 
systems as non-equilibrium dynamic systems |22]j23[1 , the chiral degree of freedom in the fully 
frustrated XY model and lattice gauge theory [j25[| . The critical exponents as well as the 
critical temperature can be extracted either from the power law behavior of the observables 
at the early times or from the finite size scaling. Compared with the non-local cluster 
algorithms, the dynamic approach does study the dynamic properties of the original local 
dynamics. 

However, it is not clear whether the short-time dynamic approach can systematically go 
beyond critical systems with a second order phase transition, even though first step approach 
or attempt has been made to the KT transitions and the spin glass transitions pp|-f29"||. For 
systems with a KT transition, for example, owing to the absence of symmetry breaking and 
to the fact that the system remains critical below Tkt, a clear signal as for a second order 



phase transition pl|Jr9|j2^| does not exists for the transition temperature Tkt- The exponent 
v and Tkt have not been determined. Standard techniques developed for the second order 
phase transitions do not apply here. On the other hand, besides the exponents, whether 
one can obtain other critical properties of the equilibrium state as especially the spatial 
correlation length from the short-time dynamics remains unknown, even for second order 
phase transitions. 

In this communication we propose a short-time dynamic approach to the KT transition 
taking the two-dimensional XY model as an example. From the short-time dynamic scaling, 
we extract the spatial correlation length of the equilibrium state. From the spatial correlation 
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length we estimate the transition temperature Tkt and the static exponent v. With Tkt at 
hand, the static exponent 77 and dynamic exponent z are obtained from power law behavior 
of the magnetization and Binder cumulant. There are no principle reasons to choose the XY 
model but only because there exist the most data of Monte Carlo simulations in equilibrium 
for comparison. 

The XY model in two dimensions is defined by the Hamiltonian 

H = ^H 3 ■ , (i) 

<ij> 

where Si = (Si jX ,Si jy ) is a planar unit vector at site i and the sum is over the nearest 
neighbors. In our notation, the coupling constant is already absorbed in the temperature. 
Large scale Monte Carlo simulations in equilibrium have been performed to understand the 
properties of the phase transition @-|5[]. The spatial correlation length £ and susceptibility x 
have been calculated in a temperature interval [0.98, 1.43] with lattice sizes up to 512. The 
results support a KT singularity for the spatial correlation length 

£(r) ~ exp (b t~ v ) (2) 

and for the susceptibility x{ T ) ~ £> 2 ~ v { T )-> with r ~ (T — T KT )/T KT being the reduced 
temperature. However, unconstrained four-parameter fits to the data do not yield completely 
satisfactory results || . The measured values of v and T KT from the data of £ and x are 
not very consistent and stable. 77 estimated from x( T ) ~ £ 2 ^( r ) is above 0.7 and too big 
compared with the theoretical prediction 77 = 0.25. The temperatures for the available data 
of £ and x are still far from the transition temperature T KT estimated to be around 0.89 to 
0.90 (for details, see Table | and Ref. @). However, simulations in equilibrium with lower 
temperatures are very difficult. 

We will demonstrate that from the short-time dynamic scaling, the spatial correlation 
length £ (r) of the equilibrium state can be extracted with relatively small lattices. This is 
because the non- equilibrium spatial length £(t, t) is small in the short-time regime of the 
dynamic evolution. Therefore, simulations can be performed at lower temperatures. 

In this paper we consider only the dynamics of model A, which is relaxational without 
energy and magnetization conservation. Starting from an ordered initial state, e.g. all 
Si = (Si iX ,Si t y) = (1,0), the system is updated at the temperature above T KT with the 
standard Metropolis algorithm. We stop updating at a certain Monte Carlo time t m and 
repeat the procedure. Total samples for average is from 800 to 1 200 for lattice size L = 256 
and above 400 for L = 512. The lattice size L = 256 is used in simulations for temperatures 
from T = 1.07 down to 0.975, while L = 512 from T = 0.97 to 0.94. Extra simulations with 
other lattice sizes confirm that our data have no visible finite size effect. 

The observable we measure is the magnetization defined as 

M (t) = 77 < E su*) > • ( 3 ) 

i 

From a general physical view point of the renormalization group transformation, the mag- 
netization M(t) is subject to a scaling form 

M(t,£(T))=t-*>' 2 *M(l,r 1 '*S(T)). (4) 
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When the temperature is at Tkt (or below), i.e. r = 0, £(r) — > oo and M(t) undergoes a 
power law decay M(t) ~ t^ v ^ 2z . However, for r > 0, the power law behavior is modified by 
the scaling function M(l,i -1 '*£(T)). This fact can be used for the determination £ z (t) and 
the exponent f]/z. 

In Fig. [I], the time evolution of the magnetization is displayed in log-log scale for different 
temperatures. We perform Monte Carlo simulations up to a time t m where there is visible 
deviation from the power law behavior. Actually, in the short-time regime of the dynamic 
evolution, the magnetization itself is more or less self-averaged. We may increase the lattice 
size without too much extra fluctuation. What restricts our simulations to very low temper- 
atures is only t m . Our longest updating time is t m = 50 000 at the temperature T = 0.94. 
To obtain the curve for T = 0.98 one needs 8 days in ALPHAstation 500 (400MHz), while 
ten times more for T = 0.94. 

Now from a scaling collapse of two curves with a pair of temperatures (Ti, T 2 ), we estimate 
the ratio £f/£f and rj/z. Here £i and £ 2 are the values of £(r) at the temperatures T\ and 
T 2 respectively. In Fig. 0, such a scaling plot is displayed for (Ti,T 2 ) = (0.955,0.965). We 
multiply the magnetization M(ti, £i) by an overall factor b a and rescale t\ to t\/b. According 
to the scaling form ([|) this rescaled M(ti,£i) is equal to M(t2, £2) if and only if b = £f/£f 
and a = r)/2z. Therefore, searching for the best fit between M(t 2 ,£ 2 ) and the rescaled 
M(t 1 ,( i i) we determine £*/£f an d T]/2z. In the figure, the circles represents the rescaled 
M(ti,ii) best fitted to M(t 2 ,£ 2 ). 

Theoretically the exponents r\ and 2; are defined at the transition temperature Tkt- The 
exponent v and the exponential singularity in Eq. (||D are defined for temperatures above 
but in the close neighborhood of Tkt- If the temperature is fairly above Tkt, in principle, 
all the exponents and the parameters b may have some dependence on the temperature. 
Usually this dependence on the temperature is neglected, otherwise the situation becomes 
too complicated ||. In our dynamic approach, we perform the scaling collapse of the mag- 
netization with two temperatures which are not too far away each other and therefore the 
dependence of rj/z on the temperature are actually considered. However, we assume v and 
bz be independent of the temperature. 

In Table [TI], the measured ratio of £i/£ 2 and exponent 7]/2z for different pairs of tem- 
peratures (T 1; T 2 ) are listed. Errors are estimated by dividing the total samples for the 
time-dependent magnetization into two groups only. The lowest temperature we reach is 
T = 0.94. For comparison, in Table |TJ available values of £f/£f from Ref. || (denoted 
by f) are also given. They are slightly smaller than our results. We observe that as the 
temperature decreases, these values of £*/£f do not increase sufficiently smoothly. Real 
errors of these data (and also our data) might be somehow bigger than given in the table. 
Taking z = 1.96 in Table | as input, the resulting r\ from r]/2z in Table |TJ is around 0.26 to 
0.31. As expected ||, we see a tendency that r\ will become around 0.25 as the temperature 
approaches Tkt- 

Now we fit the data in Table || to the exponential form in Eq. (0) and estimate the 
transition temperature Tkt, the exponent v and parameter bz. The best results are given in 
Table | in comparison with the those from simulations in equilibrium. Our results are fitted 
from a relatively lower temperature interval [0.94, 1.07], and agree well with those obtained 
in a temperature interval [0.98,1.43] in Ref. 0, in both cases of a unconstrained fit and 
a fit with a fixed v = 0.5. However, as pointed out by the authors of Refs. ||f3|, for the 
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unconstrained fit the minimum in the parameter space is not very stable in the directions 
of v and bz. It is somehow by chance that our value of v = 0.48 is so close to v — 0.47 
obtained in Ref. 0. When we vary T 2 in the fitting interval [0.94, T 2 ] from 1.07 to smaller 
values, the exponent v first drops down and then rises again after around T 2 = 1.00. These 
fluctuations very probably come from the fact that we do not have sufficient data points and 
accuracy for each data point. The estimate of Tkt is relatively stable. But the impression 
is that the value of Tkt might be slightly bigger than Tkt = 0.8942 given in the table, if 
the fitting can confidently be performed at the really close neighbourhood of Tkt- 

With the transition temperature Tkt at hand, we proceed to measure the magnetization 
M and its second moment at Tkt- To seek for the dynamic exponent z, we construct 
a Binder cumulant U = M^/M 2 — 1. Finite size scaling analysis leads to the short-time 



behavior [E6l|T9|,E 



17(f) - t d/z . (5) 

Finally an accurate value of r]/2z can be obtained from the power law decay of the magne- 
tization M(f) ~ t-^/ 2z (see Eq. (§). In Fig. |, M(f) and 17(f) have been plotted in log-log 
scale. Since our measurements only extend to t = 750, a lattice size L = 64 is sufficient. 
Total samples for average is 12 000. From the slopes of the curves in the figure we measure 
d/z and rj/2z, then calculate z and 77. The results are included in Table [| Our r\ = 0.238(4) 
coincides with the best estimate 77 = 0.235(5) in equilibrium. 

In conclusions, a dynamic approach is proposed to tackle numerically the Kosterlitz- 
Thouless phase transition. We demonstrate for the first time that not only the critical 
exponents but also the spatial correlation length of the equilibrium state can be obtained 
from the short-time dynamics. Taking the two-dimensional XY model as an example, the 
exponential divergence of the spatial correlation length is extracted from the short-time 
dynamic scaling. The transition temperature Tkt, the static exponents v and 77 as well as the 
dynamic exponent z are then estimated. Since the measurements are carried out in the short- 
time regime of the dynamic evolution, where the non- equilibrium spatial correlation length is 
small, we do not encounter difficulties of generating independent configurations. Compared 
with the simulations in equilibrium, we can perform simulations at the temperatures closer 
to Tkt- This method can in principle be applied or generalized to other kinds of phase 
transitions as second order phase transitions and spin glass transitions. 

Acknowledgements: Work supported in part by the Deutsche Forschungsgemeinschaft; 
Schu 95/9-1 and SFB 418. 
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TABLES 





rn Q/i i n7l 
[U.»4, i.u t j 


rn Q/i i n7l 
[U.y4, i.u / j 


rn qq i a qi t 


rn qq i /i qi t 


Tkt 


0.8942 


0.8926 


0.8953 


0.8914 


bz 


4.12 


3.82 


3.67 


3.38 


V 


0.48 


0.5 


0.47 


0.5 


1 KT 






0.8871 


0.8961 


* 






0.57 


0.5 


Z 


1.96(3) 








77 


.238(4) 




.235(5) 





TABLE I. Our results of the exponents and Tkt obtained in temperature interval [Ti,T2] in 
comparison with those in Ref. || (denoted by f). T KT and v* of Ref. || are from data of the 
susceptibility. The second and fourth column are the results with a fixed v = 0.5 as input. Our 
values of z and rj are measured at Tkt = 0.894 and r\ of Ref. H is estimated with finite size scaling 
and Monte Carlo renormalization group methods at also Tkt = 0.894. bz for Ref. || is calculated 
by taking our z as input. 



m,r 2 ) 


(.940, .950) 


(.950, 0.955) 


(.950, 0.960) 


(.955, .960) 


(.955, .965) 


(.960, .965) 


mi 


5.75(38) 


2.04(06) 


3.61(07) 


1.81(09) 


3.05(06) 


1.65(04) 


r]/2z 


.0680(11) 


.0695(16) 


.0690(17) 


.0682(21) 


.0674(12) 


.0650(21) 


(Ti,T 2 ) 


(.960, .970) 


(.965, .970) 


(.965, .975) 


(.970, .975) 


(0.97, 0.98) 


(0.97, 0.99) 




2.65(07) 


1.56(03) 


2.35(02) 


1.51(01) 


2.165(45) 


4.35(11) 


7]/2z 


0.0671(21) 


.0691(18) 


.0676(12) 


.0685(08) 


.0682(17) 


.0699(22) 


(Ti,T 2 ) 


(0.975, 0.98) 


(0.98, 0.99) 


(0.98, 1.00) 


(0.99, 1.00) 


(0.99, 1.01) 


(1.00, 1.01) 


ef/ei 


1.47(03) 


1.965(43) 


3.67(07) 


1.840(40) 


3.22(06) 


1.710(22) 


mi f 




1.839(84) 




1.605(61) 




1.600(39) 


7]/2z 


.0668(20) 


.0727(28) 


.0768(15) 


.0753(22) 


.0774(14) 


.0762(26) 


(Ti,T 2 ) 


(1.00, 1.02) 


(1.01, 1.02) 


(1.01, 1.03) 


(1.02, 1.03) 


(1.02, 1.04) 


(1.03, 1.04) 




2.690(29) 


1.564(08) 


2.334(21) 


1.474(10) 


2.131(20) 


1.424(07) 


t\ie 2 ] 




1.453(36) 




1.434(34) 




1.351(38) 


rj/2z 


.0784(24) 


.0781(11) 


.0768(19) 


.0773(13) 


.0777(11) 


.0754(16) 


(Ti,T 2 ) 


(1.03, 1.05) 


(1.04, 1.05) 


(1.04, 1.06) 


(1.05, 1.06) 


(1.05, 1.07) 


(1.06, 1.07) 


mi 


1.964(17) 


1.380(07) 


1.832(19) 


1.329(10) 


1.726(13) 


1.312(6) 


r]/2z 


.0778(18) 


.0835(42) 


.0777(12) 


.0705(46) 


.0781(19) 


.0900(24) 



TABLE II. The measured ratio £f/£f an< i exponent rj/2z for different pairs of temperatures. 
Values of are calculated from data in the table VIII of Ref. ||. 



7 



FIGURES 




FIG. 1. Time evolution of the magnetization in log-log scale. The temperatures are 0.94, 0.95, 
0.955, 0.96, 0.965, 0.97, 0.975, 0.98, 0.99, 1.00, 1.01, 1.02, 1.03, 1.04, 1.05, 1.06, 1.07 (from above). 



0.3 _ 
100 




FIG. 2. The scaling plot of the magnetization with a pair of temperatures (Ti,T2). The upper 
and lower solid lines correspond to temperatures T\ = 0.955 and T2 = 0.965. The circles are also 
the magnetization with T\ = 0.955 but rescaled to have the best fit with that of T2 = 0.965. 



FIG. 3. Time evolution of the magnetization M and Binder cumulant U at Tkt = 0.894 in 
log-log scale. To plot the figure, U has been multiplied by a constant 200. 
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